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ABSTRACT 

We present optical {UBVRI) and near-infrared (YJHK) photometry of the 
normal Type la supernova 2004S. We also present eight optical spectra and one 
near-IR spectrum of SN 2004S. The light curves and spectra are nearly identical 
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to those of SN 2001el. This is the first time we have seen optical and IR hght 
curves of two Type la supernovae match so closely. Within the one parameter 
family of light curves for normal Type la supernovae, that two objects should 
have such similar light curves implies that they had identical intrinsic colors 
and produced similar amounts of ^^Ni. From the similarities of the light curve 
shapes we obtain a set of extinctions as a function of wavelength which allows a 
simultaneous solution for the distance modulus difference of the two objects, the 
difference of the host galaxy extinctions, and Ry. Since SN 2001el had roughly 
an order of magnitude more host galaxy extinction than SN 2004S, the value 
of Ry = 2.15l!lo;22 pertains primarily to dust in the host galaxy of SN 2001el. 
We have also shown via Monte Carlo simulations that adding rest frame J-band 
photometry to the complement of BVRI photometry of Type la SNe decreases 
the uncertainty in the distance modulus by a factor of 2.7. A combination of 
rest frame optical and near-lR photometry clearly gives more accurate distances 
than using rest frame optical photometry alone. 

Subject headings: supernovae: individual (SN 2004S) — supernovae: individual 
(SN 2001el) — techniques: photometric — extinction: interstellar 

1. Introduction 

A Type la supernova (SN, plural SNe) is widely believed to be a carbon-oxygen white 
dwarf which explodes when its mass reaches the Chandrasekhar limit owing to mass transfer 
from a nearby companion (Livio 2000, and references therein). The explosion produces 
several tenths of a solar mass of radioactive ^^Ni, which decays to ^^Co, and finally to stable 
iron (Stritzinger et al. 2006). Since the progenitors of Type la SNe have approximately 
the same mass {Mch ~ 1.4 Mq), the explosions have approximately the same brightness at 
maximum light. Phillips (1993) first demonstrated unambiguously that the peak brightness 
of a Type la SN is related to the decline rate of the light curve. Since the publication 
of Phillips' classic paper, the common measure of the decline rate has been the number 
of 5-band magnitudes that a Type la SN fades in the first 15 days after maximum hght 
(designated Ami5(5)). A typical value is Ami5{B)— 1.10 mag. Values range from ~0.7 for 
SN 2001ay (Nugent et al. 2006, in preparation) to 1.93 for SN 1991bg (Fihppenko et al. 
1992; Leibundgut et al. 1993; PhiUips et al. 1999). 

In the past 13 years supernova researchers have noted various patterns in the light 
curves of Type la SNe. This has allowed the determination of the peak magnitudes, host 
galaxy extinction, and distances to these objects even if the data of a particular object 
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are irregularly spaced in time and do not necessarily overlap the light curve maxima. The 

three most widely used systems for fitting optical light curves are the Ami^{B) system itself 
(Hamuy et al. 1996; Phillips et al. 1999; Germany et al. 2004; Prieto, Rest, & Suntzeff 2006), 
the "stretch method" (Perlmutter et al. 1997; Goldhaber et al. 2001; Knop et al. 2003), and 
the multi-color light curve shape (MLCS) method (Riess, Press & Kirshner 1996a; Riess et 
al. 1998; Jha, Riess, & Kirshner 2006b). There is also the hybrid magnitude-color system of 
Wang et al. (2003a). 

Host galaxy extinction of Type la SNc is measured empirically by first determining 
reddening-free color loci, determining color excesses, and using empirically determined scale 
factors to convert color excesses to Ay or extinction in some other band. Lira (1995) found, 
from a small sample of Type la SNe unreddened in their hosts, that B — V colors are uniform 
and evolve in the same way from 30 to 90 days after the time of \/-band maximum. Nobili 
et al. (2003) investigated the optical colors of Type la SNe and found that all the colors 
except V — I are consistent with zero intrinsic dispersion at 35 days after T^B^ax). Wang et 
al. (2005) found that the B — V colors at 12 days after T{Bmax) provided the best method 
of determining the B — V color excess. With MLCS2k2 Jha, Riess, & Kirshner (2006b) 
adopt an intrinsic color oi B — V = 1.054 ± 0.049 at t = 35 days and assume an exponential 
distribution of host galaxy reddening with scale length E{B — V) = 0.138 mag. They also 
adopt a Bayesian prior on the distribution of Ry, with a minimum value of 1.8. Then they 
marginalize over the optical photometry to obtain the best values of Ay and Ry. 

The prime problem with using optical photometry to determine Ay along the line of 
sight to a Type la SN can be illustrated as follows. For normal dust in our Galaxy Ry = 
Ay / E{B -V) ^ 3.1 (CardeUi, Clayton, & Mathis 1989). A small uncertainty in theB-V 
color excess translates into a much larger uncertainty in Ay. The scale factor Ry is not 
necessarily equal to 3.1. The mean value most appropriate to galaxies that host Type la 
SNe is somewhere in the range 2.3 to 2.65 (Riess, Press, & Kirshner 1996b; Altavilla et al 
2004; Reindl et al. 2005; Wang et al. 2006; Jha, Riess, & Kirshner 2006b). Some very low 
values have been measured. The Type la SN 1999cl was reddened by dust with Ry = 1.55 
± 0.08 (Krisciunas et al. 2006). The dust in the host of SN 2003cg had Ry = 1.80 ± 0.19 
(Elias-Rosa et al. 2006). Our observations of SN 2006X are consistent with Ry pa 1.56.^^ 
Perhaps these very low values of Ry associated with highly reddened objects are the result 
of an interaction of the SN ejecta and the surrounding circumstellar material (Wang 2005). 

To account for non-standard reddening laws and to obtain the most robust value of the 
extinction towards a Type la SN, it is best to combine optical and IR photometry. Elias et al. 



See http://www.nd.edu/ kkrisciu/sn2006X.htnil. 
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(1985) first suggested that the V — K colors of Type la SNe might be uniform and therefore 
very useful for determining extinction for these objects. Using the standard reddening law 
of CardelU, Clayton, & Mathis (1989), Ay ^ cr Yj{V - K), where cr = 1.13 for Galactic 
(Ry = 3.1) dust, and cr = 1.07 for the extreme case of Ry = 1.55 for SN 1999cl. Thus, if 
we can determine a V — K color excess, it needs to be scaled only slightly to give us Ay. 
Using data of 8 Type la SNe, we determined zero-reddening V — J, V — H, and V — K 
loci of Type la SNe with mid-range decline rates (Krisciunas et al. 2000). Subsequently, we 
determined that the corresponding loci of Type la SNe with slow decline rates (Ami5(5) < 
1.0) are roughly 0.24 mag bluer (Krisciunas et al. 2004b). We note that V — J colors of Type 
la SNe are the least uniform of the three color indices. Also, our previous loci are valid for 
the time range —9<t< +27 days, where t is measured with respect to the time of B-band 
maximum. 

Since the Calan-Tololo SN survey in the early 1990's, extensive efforts have been devoted 
to discovering and following the light curves of Type la SNe, most notably the SN search 
carried out with the Katzman Automatic Imaging Telescope at Lick Observatory (Filippcnko 
et al. 2001). Because of the time required to observe the SNe, obtain image subtraction 
templates for those objects that require templates, and to reduce the data for many objects, 
one recent trend has been to publish papers on individual objects which are unusual, such 
as SNe 2000cx (Li et al. 2001; Candia et al. 2003) and 2002cx (Li et al. 2003). The analogy 
is that pathologists can better understand how normal human organs work by studying 
organ malfunction and disease. In this paper we discuss the astronomical case of SN 2004S, 
a normal Type la SN, the light curves of which are essentially clones of another object, 
SN 2001el (Krisciunas et al. 2003). 

SN 2004S was discovered using images taken on 2004 February 3.54 and 4.56 UT by 
Martin (2004). The supernova was located at RA = 06'^ 45"^ 43f5, DEC =-31° 13' 52'/5 
(J2000), some 47'f2 west and 2'.'5 south of the core of the Sc galaxy MCG -05-16-021 (Biggs 
2004). SN 2004S was confirmed to be a Type la supernova by Suntzeff et al. (2004) from a 
spectrum taken on 2004 Feb 6.1 UT with the CTIO 1.5-m telescope. In Fig. 1 we see the 
Suntzeff et al. (2004) spectrum, along with a spectrum of SN 2001el obtained by Wang et 
al. (2003b). Given the similarities of the optical and IR light curves and the spectra of these 
two objects, we feel fully justified in calling these two objects clones. 

According to Theureau et al. (2005), the heliocentric radial velocity of MCG -05-16- 
021 is 2808 km s~^. In the reference frame of the Cosmic Microwave Background (CMB) 
radiation the recession velocity is 2957 km s~^. Assuming a Hubble constant of 72 km s~^ 
Mpc"^ (Freedman et al. 2001), the distance modulus of SN 2004S and its host is m - M = 
33.07 ± 0.21 mag, where the uncertainty corresponds to an assumed random velocity of 300 
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km s ^. 

The Galactic reddening along the hne of sight to SN 2004S is E{B -V) ^ 0.101 ± 0.010 
mag (Schlegel, Finkbeiner & Davis 1998). Assuming Ry — 3.1 for Galactic dust, there must 
be at least 0.313 ± 0.031 mag of l/-band extinction toward SN 2004S. Given the offset of the 
SN from the core of its host (see Fig. 2), we expected the SN to be essentially unreddened 
by the host. An analysis of SN 2004S may thus provide useful information on the intrinsic 
colors of normal Type la SNe of mid-range decline rates. 

In this paper we present optical and near-IR photometry of SN 2004S. Other optical 
photometry of the object has been published by Misra et al. (2005). These authors obtained 
a light curve solution using their own measurements plus our optical data from the CTIO 
1.3-m published here. Important conclusions can be made about SN 2004S by using IR data 
as well. That the optical and IR light curves of SNe 200 lei and 2004S are so similar sheds 
light on the uniformity of explosion models of normal Type la SNe and also allows us to 
investigate the interstellar dust content of the hosts of SNe 2001el and 2004S beyond what 
is possible from a study of either object on its own. 

The paper is organized as follows. In §2 we present optical and infrared photometry 
of SN 2004S. In §3 we present optical and IR spectra of SN 2004S and compare them to 
several other Type la SNe. In §4 we consider the extinction suffered by SNe 2004S and 
2001el as two separate problems and also as a joint exercise. Using information about the 
relative uncertainties of the filter by filter photometry from §4, in §5 we carry out Monte 
Carlo simulations of extinction and reddening for a general model of the effect of the dust 
on the uncertainty of distance measurements of Type la SNe. This is relevant for future 
ground-based and space-based SN surveys. Our conclusions are presented in §6. 

2. Optical and Infrared Photometry 

Most of our photometry was taken with the CTIO 1.3-m telescope and the optical/IR 
imager ANDICAM. ANDICAM contains standard Johnson UBV filters, Kron-Cousins R 
and / filters and standard Caltech/CTIO JHK filters. Read out in 2 x 2 binning mode, 
ANDICAM gives a plate scale on the 1.3-m telescope of 0'.'369 px^^ for optical imaging and 
^.'274 px~^ for IR imaging. The optical field of view was 6^3 by 6^3, while the IR field of 
view was 2^34 by 2'M. 

ANDICAM also contains a 1.03 /xm filter (known as Y). This new filter exploits a 
previously unused atmospheric window. For more information see Hillenbrand et al. (2002). 
However, there are problems with their calibration of standard stars. From synthetic pho- 
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tometry of Kurucz model spectra spanning a range of temperatures Hamuy et al. (2006, 
Appendix C) obtained a relationship between Y — Kg colors and the published J — Kg colors 
of Persson et al. (1998) standards. The resulting F-band magnitudes for standard stars are 
typically accurate to better than 0.01 mag. 

This paper also contains some optical photometry obtained with the Siding Spring 2.3- 

m telescope and the ESO 2.2-m telescope. One night of JHK photometry was obtained 
with the 3.58-m Telescopio Nazionale Galileo (TNG) at La Palma. Two epochs of latc-time 
near-IR photometry were obtained with the number 1 Very Large Telescope (VLT) at Cerro 
Paranal, using ISAAC. 

In Fig. 2 we show an optical finding chart of the field of SN 2004S. We calibrated 
the UBVRI magnitudes of some of the fields stars on 12 photometric nights using Landolt 
(1992) standards, in particular the Rubin 149 field. The optical field star magnitudes are 
to be found in Table 1. The near-IR magnitudes of some of these stars are to be found 
in Table 2. The IR cahbration was carried out on 8 photometric nights with good seeing 
using the star P9123 from the hst of Persson et al. (1998) standards. For P9123 we estimate 
Y = 11.077 ± 0.010 from the fifth order polynomial of Hamuy et al. (2006). Our JHK 
photometry of the field stars is in very good agreement (often ± 0.01 mag) with values from 
the Two Micron All-Sky Survey (2MASS). 

We present UBVRI photometry of SN 2004S in Table 3 and the YJHK photometry 
in Table 4. The brightness and location of the SN did not require the use of host galaxy 
subtraction templates. Our CTIO 1.3-m photometry is based on aperture photometry using 
a typical aperture of radius 10 px. On nights of bad seeing a larger software aperture was 
used. 

The 12 dates marked with an asterisk in Table 3 were photometric nights at CTIO. On 
those occasions we tied the SN 2004S photometry directly to the Landolt system using the 
Ru 149 field. Given that the SN 2004S field and Ru 149 were always observed at an airmass 
difference of less than 0.2, and given that the range of color of the Ru 149 standards was 
reasonably large (—0.129 < B — V < 1.115) we feel that the calibration of the field stars 
near SN 2004S must be quite robust. On other occasions we determined the photometric 
zeropoints using the field star photometry hsted in Table 1. The color terms were determined 
from observations of the Landolt (1992) standards on photometric nights. 

The UBVRI light curves of SN 2004S arc shown in Fig. 3, and the YJHK light curves of 
SN 2004S are shown in Fig. 4. The light curve templates shown in these figures are derived 
from the photometry of SN 2001el (Krisciunas et al. 2003). We simply adjusted the hght 
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curve templates in the Y-direction to minimize the statistic of the fits to SN 20048.^'' 

We note that because the optical and IR photometry of SN 2001el was obtained with 
different telescopes and different filters, the BV photometry of Krisciunas et al. (2003) 
required correction to a uniform photometric system. We adopted the filter system of Bessell 
(1979, 1990). While the "method of S-corrections" is not perfect (Stritzinger et al. 2002), it 

resolves the systematic differences of ANDICAM B and V photometry compared to SN data 
obtained with the CTIO 0.9-m telescope. Graphical and tabular values of the S-corrections 
can be found in Krisciunas et al. (2003) and Krisciunas et al. (2004c). Officially, each SN 
should have its own set of S-corrections. While many Type la SNe have similar spectra, 
they are not identical, so applying the S-corrections adds scatter to the data. However, since 
the S-corrections we used were based on SNe 1999ee and 2001el, and SN 2001el is so similar 
photometrically to SN 2004S, correcting SN 2004S data in the same manner as we corrected 
data of SN 2001el may add only ± 0.01 mag of scatter. 

Similary, we calculated IR S-corrections to place the ANDICAM JHK photometry on 
the system of Persson et al. (1998). We found that corrections to the R- and /-band photom- 
etry actually made the photometry of SN 2001el from different telescopes more discordant, 
not less. Here we use only the ANDICAM photometry of SN 2001el (without S-corrections) 
to produce R- and /-band templates for a direct comparison of photometry of SNe 2001el 
and 2004S.^^ The key point is that we have made the same corrections to the photometry of 
the two SNe in the same way, filter by filter, so that a direct comparsion of the light curves 
of the two objects is consistent. 

The one exception is the ^-band. On the one hand, the shapes of the [/-band light ciuvcs 
of SNe 2001el and 2004S are quite similar. On the other hand, the [/-band photometry of 
SN 2004S is "too faint", on average, compared to SN 2001el. We can use a single function 
to fit the BVRIJHK data sets (see Fig. 14 below), but the mean difference of the [/-band 
data is discrepant by 0.26 ± 0.07 mag. Jha et al. (2006a, Fig. 9) have shown that normal 
Type la SNe at maximum light have B — V ^ —0.1, but that the U — B colors of these 
objects range more than half a magnitude. They emphasize that "this is not an artifact 



^^For BVRI we used cubic splines to fit the SN 2001el light curves. For UJHK we used polynomials. 

The J7-band photometry of SN 2001el from 29 to 65 days after the time of S-band maximum was synthetic 
photometry based on HST spectra kindly made available by P. Nugent. 

i^In 2001 ANDICAM was used with the CTIO 1.0-m Yale-AURA-Lisbon-Ohio (YALO) telescope, while 
beginning in 2003 it has been used on the 1.3-m ex-2MASS telescope. We assume that the optics of the 
primary and secondary mirrors of the two telescopes are sufficiently similar that no serious photometric 
differences result. We also note; that all of the [/-band photometry of SN 2001el was obtained with the CTIO 
0.9-m and 1.5-m telescopes, not with ANDICAM. 
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of the reddening correction, nor can it be explained by variation in the extinction law in 
these external galaxies." Their data imply that "there is a significant intrinsic dispersion in 
[/-band peak brightness even after accounting for variations in light-curve shape." Since the 
ultraviolet portion of the spectra of Type la SNe contains many metallic lines, a significant 
fraction of the dispersion of [/-band photometry from object to object must be due to the 
metallicity of the progenitors (Hoeflich, Wheeler, & Thielemann 1998; Podsiadlowski et al. 
2006). Suffice it to say that a full understanding of the [/-band photometry and deductions 
that may be made from it are beyond the scope of the present paper. 

A visual inspection of Figs. 3 and 4 demonstrates the remarkable similarity of the light 
curves of SNe 200 lei and 2004S. Any differences are more obvious if we plot the residuals 
of the fits, as shown in Figs. 5 and 6. Using SN 2001el templates, the reduced values of 
the fits to the CTIO 1.3-m data of SN 2004S are 1.24, 1.67, 1.12, 1.50, and 3.45 for [/, B, 
V, R, and /, respectively. For J, H, and K the reduced values of the fits are 1.90, 5.62, 
and 0.99, respectively. 

Very little F-band photometry of SNe has been pubhshed so far. Previously, we pub- 
lished some F-band data of SNe 1999ee and 2000bh (Krisciunas et al. 2004b) and noted that 
the second maximum of SN 2000bh was brighter than the first maximum. SN 2004S has 
shown the same behavior. Hamuy et al. (2006) indicate that many Type la SNe observed at 
Las Campanas Observatory as part of the Carnegie Supernova Project also show stronger 
second maxima in the y-band compared to the first maxima. 

The origin of the secondary maximum in the fight curves of Type la SNe has recently 
been discussed by Kasen (2006). His results derive from 1-D (spherically symmetric) mod- 
elling. The basic idea is that a mixture of 75 percent ^^Co, 24 percent iron, and 1 percent 
^^Ni, has a near-lR emissivity that sharply peaks at a temperature of 7000 K. This cor- 
responds to the composition and temperature of the ejecta of a Type la SN 40 days after 
explosion, or 20 days after 1{Bmax)- (The exact composition is not critical, however, so long 
as we are considering iron group elements.) At this time the iron group elements undergo 
a transition from a state of double to single ionization. This in turn causes a redistribution 
of energy from optically thick shorter wavelengths to optically thin near-IR wavelengths. If 
the y-band secondary maximum is brighter than the first maximum in most or all cases, 
the process which gives rise to the secondary maximum must be happening to the greatest 
degree at 1.03 jim. 

Kasen (2006) also suggests that a transition from singly ionized iron group elements to 
neutral atoms, occurring some 80 days after T{Bjnax), may give rise to a third maximum. 



Some preliminary light curves can be seen at littp://cspl.lco.cl/ cspuserl/CSP.html. 
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In Figs. 4 and 7 we see evidence for this phenomenon in the hght curves of SN 2004S, at 
the 2.8(j level in the J-band and up to 5.3a in the i^-band. There exist very few near-IR 
measurements of Type la SNe at these late epochs. We now have an extra movitation to 
observe these objects in the near-IR at late epochs — to check Kasen's prediction. 

Using the light curve fitting code of Prieto, Rest, & Suntzeff (2006), which relies on 
BVRI data, we find for SN 2004S that T(S„„^) = JD 2,453,039.87 ± 0.25, Ami.,{B)= 1.14 
± 0.01, and E{B — V)host = 0.00 ± 0.01. The derived value of the decline rate is statistically 
equal to that of SN 2001el, Ami5(5)= 1.13 ± 0.04 (Krisciunas et al. 2003).^^ The Ami5(5) 
solution confirms that SN 2004S had minimal host galaxy reddening, Ay(host) ^ 0.00 ± 
0.03. A more robust solution for the host galaxy extinction is obtained using the optical and 
IR data (see below). The Amis (5) method gives a distance modulus of m — M = 33.32 ± 
0.17 mag onan h — 0.72 scale (Preedman et al. 2001).^^ This is to be compared with m — M 
— 33.07 ± 0.21 mag from the recession velocity in the CMB frame. 

Using the SN 2001el l^-band template, we estimate that the F-band maximum of 
SN 2004S is Vmax = 14.48 ± 0.03. With Ay (tot) = 0.38 mag (see §4) and adopting a 
distance modulus of m — M = 33.07 ± 0.21 mag from the radial velocity in the CMB frame 
and Bubble's Law, we obtain My (max) = —18.97 ± 0.22 mag. The expected value for a 
Type la SN of average decline rate is My = —19.13 on an /i = 0.72 scale (Phillips et al. 
1999). Thus, SN 2004S has an absolute V magnitude at maximum within 1-a of the nominal 
value for its decline rate. 

Using the SN 2001el templates, we estimate that the S-corrected, K-corrected apparent 

magnitudes at maximiim light of SN 2004S were Jmax = 14.806. H,nax = 14.975, and K^nax — 
14.766, with uncertainties of ± 0.03 mag. The total extinction due to dust in our Galaxy and 
in the host of SN 2004S is Aj = 0.107, Ah = 0.068, and Ak = 0.043 mag, with uncertainties 
of ± 0.02 mag or less (see §4). Adopting a distance modulus oi m — M = 33.07 ± 0.21 
mag from the radial velocity in the CMB frame and Rubble's Law, we obtain IR absolute 
magnitudes at maximum light of Mj = —18.37, Mh — —18.16, and Mk — —18.35, with 
uncertainties of ± 0.22 mag. These are to be compared with the mean values of more 
than 20 objects with Ami5(S) < 1.8 of Mj = -18.61, M^^ = -18.28, and M;^ = -18.44 
(Krisciunas ct al. 2004c, Table 17). Taking the numbers at face value, SN 2004S is slightly 
underluminous in the IR, which contradicts the implication from MLCS2k2 analysis that 



I'^In the context of the MLCS2k2 system Jha, Riess, & Kirshner (2006b) give A = -0.19 ± 0.03 mag for 
SN 2001el. Using a copy of that software we obtain A = -0.24 ± 0.04 mag for SN 2004S. A < is a measure 
of the over luminosity of these SNe compared to the fiducial object of the system. 

18 We define /i = Hq / (100 km s"! Mpc"!). 
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the object is somewhat over luminous. In any case, the near-IR absolute magnitudes are still 
consistent with the notion that Type la SNe with Ami5(i?) ;^ 1.8 are standard candles in 
the near-IR (Krisciunas, Phillips, & Suntzeff 2004a), with an intrinsic dispersion of roughly 
±0.15 mag. 

For SN 2004S Misra et al. (2005) found Ay(tot) = 0.542 ± 0.167 mag. This implies host 
galaxy extinction of Ay(host) ~ 0.23 mag, with a rather large uncertainty. These authors 
adopt a distance modulus of 32.94 mag on an /i = 0.65 scale, equivalent to m — M = 32.72 
mag on an /i = 0.72 scale. They used a somewhat lower heliocentric redshift for the host 
galaxy (2730 km s~^) and corrected for the Virgoccntric infall of the Local Group, but do not 
quote a velocity in the CMB frame. One can also use the flow model of Tonry et al. (2000) 
to give an estimate of the distance to SN 2004S. Along a vector towards the direction of a 
galaxy in the sky, one determines the distance at which the flow model gives a radial velocity 
that matches the observed hehocentric velocity. For the host of SN 2004S that distance is 
36.25 Mpc on an /i = 0.784 scale, using the Theureau et al. (2005) heliocentric radial velocity 
of 2808 km s~^. On an h = 0.72 scale the equivalent distance modulus is m — M = 32.98 
mag, statistically in agreement with the value of 33.07 ± 0.21 mag from the CMB velocity 
and Hubble 's Law. 

Since SN 2001el is known as one of the few Type la SNe to exhibit measurable polar- 
ization, implying an asymmetrical explosion (Wang ct al. 2003b), SN 2004S may also have 
exhibited the same effect. See Chornock et al. (2006) for a discussion. 

A flnal topic to mention in this section is the late-time photometry given in Tables 3 

and 4. As Sollerman et al. (2004) found in the case of SN 2000cx, six months to a year after 
T(i?maa;) SN 2004S coutinucd to get fainter at optical wavelengths, but in the near-IR it 
reached a plateau and stayed constant. This conflrms that at late times the near-IR flux of 
a Type la SN becomes a larger and larger fraction of the total flux (see Fig. 13 of Sollerman 
ct al. 2004). We show in Fig. 7 our V- and H-hand photometry of SN 2004S. The S-band 
photometry from 50 to 410 days after T{Brnax) gives dm/dt = 1.50 ± 0.01 mag per 100 days. 
The late-time photometry {t > 196 d) in the V, R and / bands gives dm/dt = 1.42 ± 0.04, 
1.55 ± 0.12, and 1.40 ± 0.10 mag per 100 days, respectively. Sollerman et al. (2004) found 
decline rates of ~1.4 mag per 100 d for the optical light curves of SN 2000cx. Stritzinger & 
Sollerman (2007, in preparation) present late-time photometry of SN 2001el. 

At late times the optical minus IR color of SN 2004S became very red. V — H — 
+0.27 ± 0.10 at t = 76 d (see Figs. 12 and 13 below). At t ^ 3A2 d the V - H color was 
one magnitude redder (-1-1.27 ± 0.20), and at t = 416 d it was one magnitude redder still 
(+2.27 ± 0.20). 



- 11 - 



3. Optical and Infrared Spectra of SN 2004S 

In Table 5 we summarize some basic information about the spectra of SN 2004S pre- 
sented here. The spectra were obtained with the CTIO 1.5-m telescope, the 3.58-m Tele- 
scopio Nazionale Galileo (TNG), the 2.56-m Nordic Optical Telescope (NOT), and the 10-m 
Keck telescope, and were reduced following the standard recipe for longslit spectroscopy. 
Spectropolarimetry of SN 2004S, obtained with the Keck telescope, and a comparison to 
SN 2001el, are also given by Chornock et al. (2006). 

In Fig. 8 we show optical spectra of SN 2004S obtained between 1.7 and 42.5 days 
after T{Bmax)- In Fig- 9 we compare the optical spectra of SN 2004S at three epochs 
to several other normal Type la supernovae: SN 1996X (Salvo et al. 2001), SN 1998bu 
(Hernandez et al. 2000), SN 1999ee (Hamuy et al. 2002), SN 2001el (Wang et al. 2003b), 
and SN 2003du (Gerardy et al. 2004; Stanishev et al. 2006). Clearly, the spectral evolution 
of SN 2004S closely follows that of other normal Type la SNe. It is worth noting that the 
largest differences between the SNe are observed in the Ca II H&K and infra-red triplet lines. 
Otherwise, all the spectra are nearly identical. 

In Fig. 1 one can see how remarkably similar SNe 2001el and 2004S are at maximum 
light. Branch et al. (2006) investigated the maximum-light spectra of 24 Type la SNe. 
Their "core-normal" group of 7 objects included SNe 1998bu and 2001el. They noted the 
high- velocity Ca II absorption of SN 200 lei and also identified the three local minima in the 
spectrum between 4600 and 4900 A as being due to high- velocity Fe II. Since SN 2004S shows 
the same features, the identification undoubtedly holds for SN 2004S as well. Branch et al. 
(2006) also noted that SN 2001el had the strongest high-velocity Fe II of all "core-normal" 
SNe. In Fig. 1 one can see that high- velocity Fell features in SN 2004S are as strong as 
in SN 2001el, thus strengthening the similarity between the two SNe. Besides, the two SNe 
had equal Si II ratios, 7?.(SiII) ^ 0.3, as defined by Nugent et al. (1995). 

Despite the similarities between SNe 2001el and 2004S, there are a few spectroscopic 
differences. The depth of the Si II absorption observed at 6150 A is weaker in SN 2004S. 
Part of the blend due to Fe II just blueward of 5000 A is different in the two SNe. The 
greater dust extinction along the line of sight to SN 2001el produces some absorption due to 
the Na D lines at 5889/5896 A. SN 2001el had a shght shoulder in its Ca II absorption at 
7900 A whereas SN 2004S did not. 

Larger differences are observed in the Call infra-red triplet (Call IR3). The differences 
may be attributed to the presence of high- velocity features (HVFs) approximately centered 
at 20000 kms^^. In Fig. 10 we see a comparison of the velocity profiles of the Call IR3 in 
SNe 1999ee, 2001el, and 2004S. The line profiles are normalized to a local pseudo-continuum. 
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approximated by a straight line in a way similar to the definition of the Si II ratio (see Nugent 
et al. 1995). Around T{B^ax) the Call HVFs in SNe 2001el and 2004S have equal strength, 
but at the later epochs they are stronger in SN 2004S. At the same time the photospheric 
Call line is stronger in SN 2001el. The HVFs in SN 1999cc at the time of maximum light 
are weaker than in SNe 2001el and 2004S, but are stronger 20 days after. Note, however, 
that while SN 1999ee also had high-velocity Ca II, it was placed by Branch et al. (2006) in 
their "shallow-silicon" category; it had weaker Si II and S II features than SN 2001el. We 
note that SN 1999ee had a much slower dechne rate (Ami5(5)= 0.94 ± 0.06; Stritzinger 
et al. 2002), which, along with the shallower absorption fines of singly ionized species, is 
consistent with it having been a hotter explosion than SNe 2001el and 2004S. 

Another difference emerges when comparing the velocity inferred from the minimum 
of the Sill A6355 line. Around T(5„„^) SN 2004S had a velocity of ~ 9300 kms-\ while 
in SN 2001el it was ~ 10200 kms~^. Analyzing the spectrophotometric properties of a 
sample of Type la SNe, Benetti ct al. (2005) found that they can be statistically divided 
into three groups mainly on the behavior of the rate of decrease of the expansion velocity 
of the Sill A6355 absorption after the time of maximum light, {v) ('velocity gradient'), and 
also according to the velocity of the Sill A6355 at T{Bmax), '^rrtax(Sill). Respectively, the 
groups have been called high velocity gradient (HVG; (v) = 97±16 kms~M~\ Vmaxi^^^^) = 
12200±1100), low velocity gradient (LVG; (v) = 37±18kms-M-^ Vmax{SilT) = 10300±300), 
and FAINT ({v) = 87±20 kms-M-\ ^^^^(Sill) = 9200±600). The FAINT group includes 
SNe that are intrinsically dim, on average ~ 2 mag fainter than SNe belonging to the other 
two groups. SN 2001el had v — 31±5 kms~^d~^, and Benetti et al. (2005) placed it in the 
LVG group. For SN 2004S, however, we measure = 69 ± 5 kms~^d^^, which is in between 
the LVG group on one hand, and HVG and FAINT groups on the other. Along with the 
relatively low velocity at the time of maximum light, this suggests that SN 2004S may be an 
intermediate object in the context of the analysis of the Si II A6355 line velocity evolution. 

A single IR spectrum of SN 2004S was obtained 15 days after T{B.max) using the Near- 
Infrared Camera and Spectrograph at the TNG. An Amici prism was used as a disperser, 
providing the whole near-lR spectral range in one exposure but with very low resolving 
power of ~ 100. In Fig. 11 we compare the near-IR spectra of SNe 2004S, 1999ee, and 
1998bu; the near-IR spectra of these three objects are dominated by the same singly ionized 
species, mostly of iron group elements. Marion et al. (2003) showed that the abrupt change 
of the flux near 1.52 //m (rest wavelength ~ 1.57/xm) observed in many Type la SNe deflnes 
the transition from partial to complete silicon burning. In SN 2004S we measure a velocity 
of ~ 11000 km s~^, which is consistent with the measurements of other SNe and the model 
predictions for a normal- bright Type la SNe (Marion et al. 2003). 
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4. Considerations of reddening and extinction of SNe 2001el and 2004S 

The V — J, V — H, and V — K color curves of SN 2004S are shown in Fig. 12. Using 
our zero-reddening loci for Type la SNe with mid-range decline rates, we find total color 
excesses (i.e. due to Galactic dust and host galaxy dust) of E{V — J) = 0.324 ± 0.060, 
E(y -H)^ 0.371 ± 0.056, and E{V - K) ^ 0.316 ± 0.030 mag. We note that IheV -K 
curve is by far the best match as to shape. To scale the color excesses to estimates of Ay, the 
appropriate coefficients are 1.394 ± 0.110, 1.217 ± 0.058, and 1.130 ± 0.029, respectively. 
These coefficients differ slightly from the values in Eqs. 2, 3, and 4 of Krisciunas et al. 
(2004b). The Cardelli, Clayton, & Mathis (1989) values of A^/Ay are appropriate for 
determining extinction of stars with normal spectra. Since SNe have such different spectra 
we used the values in Table 8 of Krisciunas et al. (2006), which are based on spectra of 
SN 1999ee at maximum light. The uncertainites assigned to these coefficients correspond to 
a 20 percent uncertainty in Ax/ Ay. For SN 2004S we obtain three estimates of Ay, namely 
0.452 ± 0.090, 0.452 ± 0.071, and 0.357 ± 0.035 mag. The weighted mean is Ay (tot) = 
0.384 ± 0.030 mag, of which 0.313 ± 0.031 mag is due to dust in our Galaxy. It follows that 
Ay (host) = 0.071 ± 0.043 mag, implying a finite but small amount of host galaxy extinction. 
This is consistent with the host galaxy extinction derived using the Ami5(5) method. 

Since we have already demonstrated that the optical and IR light curves of SNe 2001el 
and 2004S are so similar, it is no suprise that the shapes of the photometric color curves 
of the two objects are also nearly identical. By adjusting a color template from SN 2001el 

to the corresponding color curve of SN 2004S we can determine the difference of the colors 
and therefore the difference of the color excesses. In Fig. 13 we show the same SN 2004S 
data as in Fig. 12 but with the color templates derived from data of SN 2001el, adjusted to 
minimize the statistic. We find that we need to shift the SN 2001el templates by AiV — J) 
= -0.144 ± 0.044, A(y -H) = -0.175 ± 0.035, and A(y - K) = -0.187 ± 0.068 mag. 

We may next correct for the effect of Galactic dust. SN 2001el has E(i? — V)Gai = 
0.014, while SN 2004S has E(fi - V)Gai = 0.101 (Schlegel, Finkbeiner & Davis 1998). Thus, 
the difference of the Galactic components of y-band extinction toward these two objects 
is (0.101 - 0.014) times Ry of 3.1, or 0.270 ± 0.027 mag. The differences of the color 
excesses due to dust in our Galaxy are: AE{V - J)goX = 0.270 / 1.394 = 0.194 ± 0.019, 
AE(y - H)Gai = 0.270 / 1.217 = 0.222 ± 0.022, and AE(y - K)Gai = 0.270 / 1.130 = 
0.239 ± 0.024 mag. We find that the host galaxy dust of NGC 1448 reddened SN 2001el 
more than the host galaxy dust of MCG -05-16-021 reddened SN 2004S by these amounts: 
AE(y - J)host = 0.194 + 0.144 = 0.338 ± 0.048, AE{V - H)host = 0.222 + 0.175 = 0.397 ± 
0.041, and AE{V - K)host = 0.239 + 0.187 = 0.426 ± 0.072 mag. To obtain the difference 
of host galaxy y-band extinction, we need to scale these values by coefficients appropriate 
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to the dust in the two host galaxies. Krisciunas et al. (2003) found Ay (tot) = 0.57 ± 0.05 
for SN 2001cl. With Ay(Gal) 0.043 towards SN 2001el, Ay(host) ^ 0.53 mag. Thus, 
SN 2001el suffered approximately an order of magnitude (0.53/0.07 ~ 7.6) more host galaxy 
extinction compared to SN 2004S. For SN 2001el Wang et al. (2003b) give Ry = 2.88 ± 0.15 
from the peak wavelength of the interstellar polarization.^^ Jha, Riess, & Kirshner (2006b) 
give Ry = 2.40 ± 0.19 for the host galaxy dust of SN 200 lei. Using the values in Table 8 
of Krisciunas et al. (2006) and this range of values of Ry, we find that the difference of host 
galaxy V-band extinction of the two SNe is Ay = 0.471 ± 0.027 mag. Adding this to the 
magnitude difference of the y-band light curves, corrected for Galactic extinction (column 
5 of Table 6 below), the imphed difference of the distance moduh of the two objects is 1.936 
± 0.044 mag. 

It is not absolutely necessary to make any a priori assumptions about the parameters 
of interest here. In Table 6 we lay out a spread sheet relating to the implied dust extinction 
in the hosts of SNe 2001el and 2004S. In column 2 we give the magnitude offsets necessary 
to minimize the statistic of fitting the SN 2001el hght curve templates to SN 2004S. In 
column 5 of Table 6 we have corrected the values of column 2 for the extinction due to 
Galactic dust. 

In Fig. 14 we show the values from column 5 of Table 6 vs. the wavelengths of the 
photometric filters. As the wavelength tends to infinity, the effect of the dust in the two host 
galaxies becomes negligible. In other words, in the absence of host galaxy dust extinction, 
all the points (excepting U) would lie along some horizontal line near the top of the graph. 
This is most directly interpreted as the difference of the distance moduli (= A/x) of the two 
objects. ^'^ 

We have carried out a multi-dimensional minimization, to find the values of Ry, Ay 
and A/i that give the best match to the values in column 5 of Table 6. See Table 7 and Fig. 
15. We find Ry = 2.lbtoll, Ay = 0.472 ± 0.025 mag, and A/x = 1.936 mag. Note the close 
agreement of Ay and Ayu with the values obtained from a consideration of the V—[JHK\ 
color curves of SNe 2001el and 2004S. 



^^In our view the polarization curve shown in Fig. 8 of Wang et al. (2003b) does not have a well constrained 
maximum, so their value of Ry should probably be assigned a larger numerical uncertainty. 

^^Because Type la SNe exhibit an intrinsic dispersion at any given decline rate of ± 0.15 mag, two objects 
with identical explosion mechanisms, or at least producing the same amount of ^^Ni, would not necessarily 
give the exact same extinction-corrected absolute magnitudes. And even if the two objects produced the 
same number of ergs per second, the explosions could be asymmetric. One elongated fireball might be 
observed end-on, while the other could be observed side-on, giving different apparent brightnesses. 
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Using only BVRI data we find Ry = 2.511^;^^. This is statistically in agreement with 
the values of Jha, Ricss, & Kirshner (2006b) and Wang ct al. (2003b) for the host galaxy 
dust of SN 200 lei. However, our value derived from BVRIJHK data is considerably smaller 
than that of Wang et al. (2003b), derived from polarization data. Our "best" value implies 
that the dust in the host of SN 2001el had an Ry value considerably less than the canonical 
Galactic value of 3.1. 

In Fig. 15 we show the 1-a contours of our extinction solutions using 4, 5, and 7 filter 
photometry. As the number of filters and the wavelength range of the filters increase, the 
contours become smaller and rounder. In Table 7 one can sec how the uncertainties decrease 
as the number of filters increases from 4 to 6. Adding the il'-band data changed the solution 
only slightly This is undoubtedly due to the larger uncertainty of the magnitude shift 
necessary to fit the SN 2004S X-band data with the corresponding SN 2001el template. 

We can estimate the total y-band extinction of SN 200 lei as follows. The extinction 
due to Galactic dust along the fine of sight is Ay (Gal) = 0.043 ± 0.004 mag. The host 
galaxy V-h&nd extinction of SN 2001el is 0.472 ± 0.025 mag (difference of SNe 2001el and 
2004S) plus the host galaxy extinction of SN 2004S (0.071 ± 0.043 mag from VJHK data, 
but possibly zero if we rely on the Ami5(B) solution using BVRI data). For SN 2001el 
Av(tot) = 0.586 ± 0.050 mag. This compares very well with our previously published value 
of 0.57 ± 0.05 mag (Krisciunas et al. 2003). 

5. Modelling distance uncertainty and dust 

For Type la SNe to be useful as accurate probes for cosmology it is important that 
systematic errors be minimized. Future probes of Dark Energy properties will need to control 
systematic uncertainties to better than a few percent to provide useful constraints (Miknaitis 
et al. 2006, in preparation). For a typical E{B — V)^ 0.10 mag, a systematic error of 0.2 in 
Ry leads to a systematic error of 0.02 mag in the distance modulus. The range of extinction 
properties seen in nearby supernovac, such as SN 2004S, suggests that it is important to 
measure the extinction law to avoid making a systematic error in distance determination. 
But our photometry of SN 2004S shows how difficult it is to precisely constrain Ry using 
just four optical filter bands. 

To better understand how well dust extinction can be corrected in a sample of well- 
observed Type la SNe, we have simulated the process used to recover Ry, Ay and A/j, for 
SN 2004S. In the simulation we have assumed the following: 

• Supernovae have a probability distribution for Ay like that found by Jha, Riess, & Kir- 
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shner (2006b). This is an exponential probability distribution function (pdf) peaking 
at zero extinction and having a scale factor of 0.4 mag. 

• Rv has a Gaussian pdf with a mean value of 2.4 and standard deviation of 0.4. 

• The extinction as a function of bandpass is given by Aj/Ay = Oj + ftj/Ry where i 
corresponds to one of the BVRIJHK bands and the coefficients and bi are given in 
Table 8 of Krisciunas et al. (2006). 

We then did a Monte Carlo simulation selecting pseduo-random values of Ry and Ay 
which define the extinctions in the BVRIJHK bands. How accurately these input parame- 
ters can be recovered depends on two factors: the number of filters used in the observations 
and the signal-to-noise ratio of the observations. The simulations were done assuming obser- 
vations in four bands (BVRI), five bands (BVRIJ), and all seven filters {BVRIJHK). The 
signal-to-noise ratio in each band was simulated by setting a variance, (Tp^^^, about the true 
extinction. Each realization of the Monte Carlo simulation selected a random value that had 
a Gaussian distribution with that variance and centered on the true bandpass extinction. 

The resulting simulated measurements were then fit with three free parameters, Ry, 
Ay, and A/i, using a minimization technique. A/i is the offset in the recovered distance 
modulus and the measure for how well the various filter set and photometric uncertainties 
perform. 

A sample result for 30,000 events with a photometric error of 0.02 mag and recovery in 
BVRIJ bands is shown in Fig. 16. The difference between the input Ay and recovered Ay 
(i.e. the Ay error) is plotted along with the corresponding Ry error to show the correlation 
between these parameters. Note that on average the recovered values of Ry and Ay have no 
significant systematic error. 

In Fig. 17 we explore how adding IR bands improves the distance estimation. As ex- 
pected, the use of BVRIJHK comes closest to the optimal recovery of the distance modulus. 
The use of only BVRI bands means that the recovered distance uncertainty rises sharply 
with photometric error. For BVRI , photometric errors of only 0.03 mag result in distance 
uncertainites of 0.2 mag, which is the same order as the intrinsic scatter in the luminosities 
of Type la SNe. We note that our simulation was quite simple and assumes that the true 
color distribution for unextincted Type la SNe at all decline rates is known with perfect 
precision. 

The addition of the rest frame J-band to the filter complement decreases the distance 
uncertainty by a factor of 2.7. The use of rest frame J-band observations greatly improves 
the precision in estimating the distance modulus when independently fitting a two-parameter 
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{Rv, Ay) extinction law. Adding H and K further improves the uncertainty of the distance 
modulus, but only by an additional 30 percent. 

Any future SN survey to be carried out with a satellite-borne telescope would wisely 
observe Type la SNe in rest frame optical and IR bands. To determine the cosmic equation of 
state parameter (w) to within 10 percent requires controlling photometric systematic errors 
at the 0.01 to 0.02 mag level (Miknaitis et al., 2006, in preparation). At minimum we must 
know the appropriate mean value of Ry for high redshift SNe. Better still would be the 
option to determine Ry accurately for a large fraction of those yet-to-be discovered objects. 

6. Conclusions 

In this paper we provided optical and IR photometry of the normal Type la SN 2004S. 
Using optical and ncar-IR data, it was found to have minimal host galaxy extinction, with 
Ay (host) = 0.07 ± 0.04 mag. 

There are remarkable similarities of the shapes of the optical and IR light curves of SNe 
2004S and 2001el. This is the first time we have seen two objects exhibit such similar JHK 
light curves at the time of maximum light until 65 days afterwards. We might assume that 
the two objects had identical explosion mechanisms, or at least produced the same amount 
of ^^Ni. Given that Type la SNe exhibit an intrinsic dispersion of roughly ±0.15 mag at any 
given decline rate, SNe 200 lei and 2004S could have produced identically shaped filter- by- 
filter light curves, yet could have had absolute magnitudes that differ by 0.15 mag or more. 
For our purposes it is reasonable to interpret the difference of their extinction-corrected 
magnitudes in any band (except U) as a measure of the difference of their distance moduli. 
More generally, we may assume that they had identical extinction-corrected colors. 

SNe 2001el and 2004S not only had similar hght curves, but very similar optical and 
near-IR spectra, exhibiting high- velocity absorption due to singly ionized Ca and Fe. Except 
for the high velocity features in their ejecta prior to t ^ -|-14 d, they can be considered 
spectroscopically normal (Branch et al. 2006). 

The similarities of the light curves of SNe 2001el and 2004S allow us to use the photome- 
try, filter by filter, to derive the difference of the host galaxy extinctions without making any a 
priori assumptions about the unreddened colors of Type la SNe, values of Ry, or the redden- 
ing law between the K-band and infinite wavelength. Since SN 2001el had roughly an order 
of magnitude more host galaxy extinction compared to SN 2004S, our multi-dimensional 
minimization primarily tells us something about the dust in the host of SN 200 lei. For this 
object we find Ry = 2.151q;22 Ay(tot) ~ 0.59 ± 0.05 mag. Because of the advantages 
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of combining optical and IR photometry to determine extinction, it is important that we 
know the shapes and zeropoints of the color curves of Type la SNe. Our analysis gives us 
unreddened solutions for SNe 2001el and 2004S. 

The V — K color curves of SNe 2001el and 2004S have the same shape as our template 
based on the not-so-well sampled hght curves of 8 mid-range dechners. Our V — J and V — H 

templates do not have quite the same shapes as the color curves of SNe 2001el and 2004S. 
In the near future it should be possible to construct more robust color curves from a larger 
database of objects, which will reflect more realistically the inherent scatter exhibited by 
Type la SNe. 

We find that SN 2004S had VJHK absolute magnitudes at maximum light within 1-a 
of the mean of other Type la SNe of comparable decline rate. Any systematic differences 
may only be due to the unknown peculiar velocity of the host galaxy. The JHK absolute 
magnitudes at maximum of SN 2004S are consistent with the notion that in the near-IR 
Type la SNe (excepting the fastest decliners) are standard candles. 

We can derive host galaxy extinction much more accurately using optical and IR pho- 
tometry than using optical photometry alone, as we have shown graphically and quanti- 
tatively. Our Monte Carlo simulations indicate that rest frame BVRIJ photometry gives 
distance uncertainties 2.7 times smaller than just using rest frame BVRI photometry. As- 
suming a fixed value of Ry for high-redshift SNe runs the risk of systematic errors if the fixed 
value is far from the true mean. Here we showed that adding one additional rest frame band 
in the near-IR substantially improves the resulting uncertainty of determining the distance 
moduli of Type la SNe. Rest frame J-band corresponds to 3.6 //m in the observer frame 
at 2; ~ 2, which is observable with space-based telescopes. Having rest frame optical and 
near-IR data for Type la SNe will lead to stricter constraints on models of Dark Energy. 



This work is partly based on observations collected at the Nordic Optical Telescope 
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observed SN 2004S. 

The CTIO 1.3-m telescope is operated by the Small and Moderate Aperture Research 
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Table 1. Optical Photometric Sequence near SN 20048*^ 



StarlDi^ U B V R I 

1 16.841 (0.011) 16.442 (0.008) 15.612 (0.006) 15.140 (0.006) 14.727 (0.005) 

2 19.282 (0.066) 18.504 (0.012) 17.441 (0.010) 16.777 (0.015) 16.222 (0.008) 

3 18.134 (0.026) 18.157 (0.014) 17.609 (0.010) 17.257 (0.011) 16.884 (0.015) 

4 19.081 (0.054) 17.897 (0.011) 16.527 (0.007) 15.642 (0.005) 14.845 (0.006) 

5 19.255 (0.051) 18.556 (0.017) 17.509 (0.009) 16.875 (0.006) 16.325 (0.007) 

6 16.377 (0.008) 16.282 (0.008) 15.590 (0.005) 15.175 (0.006) 14.771 (0.006) 

7 18.206 (0.042) 17.848 (0.010) 16.996 (0.007) 16.501 (0.008) 16.056 (0.006) 

*The numbers in parentheses are 1-a uncertainties (mean errors of the mean). 
''The identifications are the same as in Fig. 2. 



Table 2. Infrared Photometric Sequence near SN 2004S* 

Star ID'' Y J H K 

1 14.450 (0.005) 14.166 (0.009) 13.794 (0.009) 13.682 (0.017) 

2 15.828 (0.010) 15.448 (0.010) 14.879 (0.009) 14.751 (0.039) 

3 16.595 (0.019) 16.441 (0.023) 16.117 (0.024) 

4 14.340 (0.005) 13.918 (0.008) 13.273 (0.009) 13.100 (0.017) 

^The numbers in parentheses are l-cr uncertainties (mean errors of the 
mean) . 

''The identifications are the same as those in Table 1 and Fig. 2. 
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Table 3. UBVRI Photometry of SN 2004S' 



JD-2,450,000 




U 




B 




V 




R 




I 




Telescope'' 


3042.72* 


14 


645 


(0.053) 


14 


554 


(0.024) 


14 


447 


(0.024) 


14 


385 


(0.022) 


14 


762 


(0.021) 


1 


3045.75* 


14 


965 


(0.060) 


14 


723 


(0.023) 


14 


491 


(0.019) 


14 


446 


(0.025) 


14 


839 


(0.022) 


1 


3048.76 


15 


317 


(0.085) 


14 


948 


(0.029) 


14 


592 


(0.034) 


14 


596 


(0.015) 


15 


015 


(0.024) 


1 


3052.72 


15 


865 


(0.061) 


15 


335 


(0.026) 


14 


805 


(0.024) 


14 


823 


(0.023) 


15 


246 


0.023) 


1 


3056.68* 


16 


452 


(0.051) 


15 


799 


(0.025) 


15 


050 


(0.022) 


14 


962 


(0.021) 


15 


228 


(0.023) 


1 


3060.59 


16 


886 


(0.054) 


16 


231 


(0.025) 


15 


256 


(0.022) 


15 


040 


(0.032) 


15 


128 


(0.022) 


1 


3066.67* 


17 


406 


(0.061) 


16 


744 


(0.026) 


15 


546 


(0.028) 


15 


161 


(0.026) 


14 


993 


(0.022) 


1 


3070.60* 


17 


729 


(0.080) 


17 


093 


(0.031) 


15 


808 


(0.022) 


15 


415 


(0.022) 


15 


125 


(0.028) 


1 


3074.57* 


17 


881 


(0.066) 


17 


357 


(0.026) 


16 


103 


(0.025) 


15 


723 


(0.023) 


15 


394 


(0.022) 


1 


3081.64* 


18 


094 


(0.063) 


17 


600 


(0.023) 


16 


377 


(0.033) 


16 


035 


(0.022) 


15 


750 


(0.023) 


1 


3087.60 


18 


330 


(0.087) 


17 


806 


(0.019) 


16 


602 


(0.015) 


16 


306 


(0.011) 


16 


105 


(0.019) 


1 


3090.55* 


18 


445 


(0.047) 


17 


841 


(0.023) 


16 


693 


(0.022) 


16 


402 


(0.021) 


16 


247 


(0.020) 


1 


3092.94 


18 


301 


(0.104) 


17 


858 


(0.046) 


16 


806 


(0.037) 


16 


473 


(0.067) 


16 


333 


(0.053) 


2 


3096.54 


18 


528 


(0.104) 


17 


901 


(0.052) 


16 


856 


(0.050) 


16 


614 


(0.028) 


16 


452 


(0.025) 


1 


3102.55* 


18 


581 


(0.147) 


17 


964 


(0.039) 


16 


967 


(0.027) 


16 


709 


(0.027) 


16 


711 


(0.034) 


1 


3109.51 


18 


786 


(0.066) 


18 


101 


(0.026) 


17 


177 


(0.020) 


16 


967 


(0.021) 


16 


954 


(0.061) 


1 


3116.50 


18 


834 


(0.061) 


18 


162 


(0.022) 


17 


325 


(0.019) 


17 


190 


(0.014) 


17 


203 


(0.052) 


1 


3122.48* 








18 


368 


(0.066) 


17 


573 


(0.032) 


17 


407 


(0.044) 


17 


299 


(0.045) 


1 


3130.51* 


18 


790 


(0.321) 


18 


259 


(0.067) 


17 


685 


(0.033) 


17 


638 


(0.044) 


17 


667 


(0.053) 


1 


3137.47 


19 


014 


(0.228) 


18 


438 


(0.045) 


17 


811 


(0.022) 


17 


799 


(0.018) 


17 


894 


(0.040) 


1 


3151.47 








18 


619 


(0.038) 


18 


122 


(0.034) 


18 


214 


(0.098) 


18 


194 


(0.036) 


1 


3238.31 














19 


778 


(0.058) 


20 


369 


(0.090) 


19 


489 


(0.089) 


2 


3267.26 














20 


394 


(0.049) 


21 


034 


(0.080) 


20 


547 


(0.151) 


2 


3357.75 








21 


731 


(0.035) 


21 


471 


(0.041) 


22 


271 


(0.123) 


21 


306 


(0.102) 


3 


3361.14 














21 


689 


(0.094) 


22 


292 


(0.247) 


21 


217 


(0.126) 


2 


3453.51 








23 


198 


(0.140) 


23 


084 


(0.117) 














3 



^Thc nights marked by an asterisk were photometric at CTIO. On these occasions the SN 2004S photometry was 
tied directly to the Landolt (1992) system using the stars in Ru 149. The numbers in parentheses are 1-it errors (mean 
errors of the mean) and include the RMS error of the photometric calibration from the standards and the scatter of 
the multiple measurements of SN 2004S with respect to their mean. 

''I = CTIO 1.3-m. 2 = Siding Spring 2.3-m. 3 = La Silla 2.2-m. 
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Table 4. Near Infrared Photometry of SN 20048^ 



JD-2,450,000 Y J H K Telescope^ 



3042 


72 


15 


170 


(0.026) 


14.945 


(0 


021) 


15.064 


(0 


022) 


14.874 


(0.059) 


1 


3045 


75 


15 


340 


(0.026) 


15.238 


(0 


024) 


15.090 


(0 


021) 


14.976 


(0.060) 


1 


3048 


76 


15 


472 


(0.026) 


15.641 


(0 


029) 


15.228 


(0 


026) 


15.116 


(0.066) 


1 


3052 


72 


15 


499 


(0.026) 


16.156 


(0 


036) 


15.160 


(0 


021) 


15.042 


(0.063) 


1 


3053 


43 








16.124 


(0 


010) 


15.151 


(0 


012) 


14.946 


(0.019) 


2 


3056 


68 


15 


299 


(0.025) 


16.258 


(0 


040) 


15.107 


(0 


021) 


14.886 


(0.048) 


1 


3060 


59 


15 


181 


(0.023) 


16.227 


(0 


040) 


15.041 


(0 


021) 


14.957 


(0.054) 


1 


3070 


60 


14 


750 


(0.020) 


15.866 


(0 


036) 


15.050 


(0 


024) 


15.007 


(0.055) 


1 


3074 


57 


14 


782 


(0.020) 


15.932 


(0 


035) 


15.260 


(0 


027) 


15.204 


(0.060) 


1 


3087 


60 


15 


544 


(0.037) 


17.013 


(0 


093) 


16.017 


(0 


046) 


16.017 


(0.148) 


1 


3090 


55 


15 


727 


(0.031) 


17.233 


(0 


099) 


16.221 


(0 


056) 


16.065 


(0.153) 


1 


3096 


54 


16 


059 


(0.055) 


















1 


3102 


55 


16 


293 


(0.043) 


18.200 


(0 


267) 


16.779 


(0 


085) 


16.561 


(0.254) 


1 


3109 


51 








18.587 


(0 


259) 


16.851 


(0 


077) 






1 


3116 


50 








18.652 


(0 


231) 


17.030 


(0 


099) 






1 


3130 


51 














18.407 


(0 


320) 






1 


3341 


81 














20.57 


(0 


18) 






3 


3342 


80 








21.35 


(0 


06) 












3 


3414 


71 








21.27 


(0 


13) 












3 


3415 


62 














20.63 


(0 


18) 






3 



The numbers in parentheses are l-cr uncertainties (mean errors of the mean). 
1 = CTIO 1.3-m. 2 = TNG. 3 = VLT number 1. 
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Table 5. Spectra of SN 2004S 



UT Date 


Telescope^ 


JD-2,450,000 


Phase*^ 


Range (A) 


Resolution (A) 


Exptime (sec) 


FcbOe.lO 


1 


3041.60 


1.7 


3100-9650 


17 


900 


Febll.96 


2 


3047.46 


7.6 


3230-10400 


25 


900 


FeblT.OO 


3 


3052.50 


12.6 


4880-7940 


20 


2x600 


Febl7.85 


3 


3053.35 


13.5 


3260-9050 


20 


2x600 


Febl8.45 


4 


3053.95 


14.1 


3360-10380 


10-12 


2x180 


Feb22.92 


3 


3058.42 


18.5 


3190-9050 


20 


2x600 


Marl 1.88 


3 


3076.38 


36.5 


3200-9040 


20 


2x900 


Marl7.85 


3 


3082.35 


42.5 


3210-9060 


15 


2x900 



^1 = CTIO 1.5-m. 2 = 3.58-m Telescopio Nazionale Galileo. 3 = 2.56-m Nordic Optical Telescope. 
4 = Keck 10-m, using LRIS. 

''Days since T{Bmax) = JD 2,453,039.87. 
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Table 6. Adjusting SN 2001el templates to the light curves of SN 2004S 



Filter 


Amagi^ 






Amag2'^ 


U 


1.725 (0.062) 


1.5288 


0.413 (0.041) 


1.312 (0.074) 


B 


1.608 (0.025) 


1.3147 


0.355 (0.036) 


1.253 (0.044) 


V 


1.735 (0.022) 


1.0000 


0.270 (0.027) 


1.465 (0.035) 


R 


1.781 (0.019) 


0.8299 


0.224 (0.022) 


1.557 (0.029) 


I 


1.859 (0.022) 


0.6050 


0.163 (0.016) 


1.696 (0.027) 


J 


1.888 (0.033) 


0.2828 


0.076 (0.008) 


1.812 (0.034) 


H 


1.919 (0.024) 


0.1785 


0.048 (0.005) 


1.871 (0.025) 


K 


1.928 (0.063) 


0.1154 


0.031 (0.003) 


1.897 (0.063) 



'^Number of magnitudes the SN 200 lei templates must be 
shifted to produce the best fits to the SN 2004S data. 

''Prom column 4 of Table 8 of Krisciunas et al. (2006). 

"^Adjustments for Galactic extinction, equal to 0.270 mag times 
the coefficients of column 3. We assumed Ry = 3.1 for Galactic 
dust and uncertainties of 10 percent. 

Values of column 2 minus those in column 4. These are the 
equivalent light curve template shifts after correcting for the ef- 
fects of Galactic extinction. 
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Table 7. Extinction solutions^ 



Filters 


Ry 










BVRI 


2.51 (+0.50,- 


-0.36) 


0.534 (+0.075,- 


-0.057) 


1.994 


BVRIJ 


2.02 (+0.30,- 


-0.28) 


0.451 (+0.041,- 


-0.043) 


1.916 


BVRIJH 


2.14 (+0.24,- 


-0.22) 


0.470 (+0.027,- 


-0.026) 


1.934 


BVRIJHK 


2.15 (+0.24,- 


-0.22) 


0.472 (+0.025,- 


-0.025) 


1.936 



^Numbers in parentheses are upper and lower l-cr errors. 

•^Difference of F-band extinction of SNe 2001el and 2004S, in mag- 
nitudes. 

^Implied difference of distance moduli of SNe 2001el and 2004S, in 
magnitudes. 
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Fig. 1.— Maximum light spectra of SN 2001el (2001 October 1), obtained with the VLT + 
FORSl, and SN 2004S (2004 February 6), obtained with the CTIO 1.5-m. The SN 2004S 
spectrum is shifted by an artibrary amount for display purposes. 

Fig. 2.— MCG -05-16-021, SN 2004S, and the field stars nearby. This by 3'.1 finder was 
made from a V-band exposure obtained with the CTIO 1.3-m telescope on 2004 February 
17 UT. (Note: a different figure will appear in the AJ version of this paper.) 

Fig. 3.— UBVRI photometry of SN 2004S from 3 to 110 days after T(5„„^). The circles 
are data obtained with the CTIO 1.3-m telescope, while the triangles are data obtained with 
the Siding Spring 2.3-m telescope. The templates are based on the hght curves of SN 2001el 
and have been adjusted in the Y-direction to minimize the statistic of the fits to the data 
of SN 2004S obtained with the CTIO 1.3-m. 

Fig. 4. — YJHK photometry of SN 2004S to 90 days after maximum light. The circles are 
data obtained with the CTIO 1.3-m telescope, while the triangles are data from one night 
with the 3.58-m TNG. The sohd black lines are fits to the data of SN 2001el, adjusted in 
the Y-direction to minimize the statistic of the fits to the CTIO 1.3-m data. The dashed 
blue lines are estimates of a linear decline of the light curves. At i = 69 and 76 days the 
if-band data are 4.3(7 and 5.3(7 above the dashed line. At t = 76 days the J-band datum is 
2.8cr above the dashed line. This late time photometry is evidence that SN 2004S exhibited 
a third hump in its light curve, in agreement with the prediction of Kasen (2006). 

Fig. 5. — Residuals of the UBVRI data of SN 2004S and the light curve templates based on 
fits to the data of SN 2001el. "A" is in the sense "data minus templates". 

Fig. 6. — Residuals of the JHK data of SN 2004S and the light curve templates based on 
fits to the data of SN 2001el. "A" is in the sense "data minus templates". 

Fig. 7. — All of our V- and i/-band of SN 2004S from Tables 3 and 4. The photometric 
errors are smaller than or equal to the size of the points. 

Fig. 8. — Optical spectra of SN 2004S from 1.7 to 42.5 rest frame days after the time of 
maximum light. The number of days since B-band maximum is given at the right hand side. 

Fig. 9. — Comparison of optical spectra of various Type la SNe and SN 2004S. 

Fig. 10. — Evolution of the high-velocity blue shifted Ca absorption from 1 to 18 days after 
the time of maximum light of SNe 1999ee, 2001el, and 2004S. 
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Fig. 11. — Comparison of near-IR spectra of the Type la SNe 2004S, 1999ee, and 1998bu at 
15 days after maximum light. 

Fig. 12.— V - J, V - H, and V - K color curves of SN 2004S. The sohd hues are the 
zero-reddening loci based on eight Type la SNe of Krisciunas et al. (2000), adjusted by the 
corresponding color excesses to minimize x^- 

Fig. 13. — y — J, V — H, and V — K color curves of SN 2004S. The solid lines are fits 
to the corresponding data of SN 2001el and adjusted in the Y-direction to minimize the 
statistic. 

Fig. 14. — Values from column 5 of Table 6 vs. the wavelengths of the photometric bands. 
In the limit as the wavelength tends to infinity we obtain a measure of the difference of the 
distance moduli of SNe 200 lei and 2004S (the horizontal dashed line). The (red) solid line 
is the (Rv = 2.15) solution that minimizes the statistic of the fit. The (blue) dashed line 
shows that normal Galactic dust with Ry = 3.1 does not fit the data points as well. 

Fig. 15. — Contours in the Ay-Ky plane. The X-axis represents possible values of the 
difference of the host galaxy y-band extinction of SN 2001el vs. SN 2004S. The Y-axis 
represents the value of Ky pertaining primarily to the dust in the host of SN 2001el. We 
plot the 68 percent contour levels for three solutions, using 4, 5, and 7 filter photometry, 
respectively. 

Fig. 16. — Correlation of uncertainty in Ry vs. uncertainty in Ay. We have simulated 
the recovery of the host galaxy extinction for 30,000 SNe. For this simulation wc picked 
at random values of Ry and Ay, then recovered them using BVRIJ photometry. Here we 
assumed a y-band photometric accuracy of ± 0.02 mag. Other bands had correspondingly 
larger or smaller errors on the basis of our actual SN 2001el and SN 2004S data. 

Fig. 17. — The uncertainty in the distance modulus for Type la SNe as a function of the 
photometric accuracy and the number of filters used to determine the photometric solution. 
The dotted line corresponds to the photometric accuracy. It is not possible to recover a 
distance modulus more accurately than the level of accuracy of the photometry. 
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Krisciunas et al. Fig. 1 
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Krisciunas et al. Fig. 3 
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Krisciunas et al. Fig. 6 
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